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2. Prediction by Numerical Methods of Self-Sustained Oscillations 
in a Two-Axis Model of the Nonlinear Sampled-Data LST System 

2-1. Introduction 

It has been demonstrated in [2] that the stability equations as 
a result of the describing function analysis may be solved by a 
numerical -iterative technique. In Chapter 1, the numerical technique 
is applied to the prediction of self-sustained oscillations in a two- 
axis model of the nonlinear LST system with continuous-data. In this 
chapter the same numerical technique will be extended to the prediction 
of self-sustained oscillations in a two-axis model of the nonlinear 
LST system with sampled-data . 

In the two-axis model of the LST system, the system contains two 
nonlinearities, and the general form of the stability equation in the 
sampled-data case is 

1 + G B (T,n)N(A,n,<t>) + G c (T,n)N 2 (A,n,4>) = 0 (2-1) 

where Gg(T,n) and G^(T,n) are transfer functions which describe the 
linear elements of the coupled systems, and N(A,n,<j>) represents the 
discrete describing function of the CMG nonlinearity. The amplitude of 
the sinusoidal signal at the input of the nonlinearities is denoted by 
A; T is the sampling period in second, and n is the integer which relates 
the period of the sustained oscillation to the sampling period; <j> is the 
phase angle of the input sinusoid. Therefore, 



- nT = period of sustained oscillation 


( 2 - 2 ) 


or 

2n 

0J C = nT = frequency of sustained oscillation in rad/sec (2-3) 

2-2. The Sampled-Data Two-Axis L$T System Model and Its Stability 

Equation 

Figure 2-1 illustrates the signal flow graph representation of 
the sampled-data single-axis 1ST system. The nonlinear friction of 
the CMG is modeled by the branch with the gain N, where strictly, N 
denotes the discrete describing function N(A,n,<J>). The sample-and-hold 
at the input of the nonlinearity is added for analytical purposes, so 
that the overall system may be treated by the z-transform method. 

Since the actual system does not have the sample-and-hold preceding N, 
the system configuration of Fig. 2-1 should be considered as an approxi- 
mation to the sampled-data single-axis LST system. 

Figure 2-2 shows the two-axis model of the LST system, with the 
coupling made at the output stages. The coupling coefficients « 2 and 
may be varied for various degrees of coupling between the two single- 
axis systems. 

For the study of self-sustained oscillations, we set x-j = x 2 = 

From the signal flow graph of Fig. 2-2, the following cause-and-effect 
relationships are written for the variables at the outputs of the four 
samplers: 

E* = -A*E* + A*N*X* + A*E* - A*N*X* (2-4) 

E 2 = “ A T E 2 + A 3 N * X 2 + A 2 E T " A 4 N * x * 


(2-5) 



K s + K t 

V s > ■ -" 4- 1 

G-, (s) = K q + KjS 

^6^ = J~ s 

K I 

G 2^ K q s + Kj 

:J 7 (s) ' TT*' 

V 

1 n- TS 

G. (s) = 1 " e 
ho v ; s 


Figure 2-1. The single-axis 

sampled-data LST system model. 
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XJ = -B*N*X* + B*E* 
X* = -B*N*X* + B*E| 



( 2 - 6 ) 

(2-7) 

( 2 - 8 ) 

(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 

( 2 - 12 ) 

(2-13) 

(2-14) 

(2-15) 

(2-16) 


In the last nine equations all functions are functions of s. The 
sampled signal flow graph of the system isdrawn in Fig. 2-3 using 
Eqs. (2-4) through (2-7). 

The determinant of the sampled signal flow graph is obtained from 
Fig. 2-3 by use of Mason's gain formula. The expression is given in 
terms of the z-transform functions: 
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A(z) = [1 + 2A-j (z) - a|(z)’ + A^(z)] + 2N{A,n,<J))[B 1 (z) - B 2 (z)A 3 (z) 

+ A 2 (z)B 2 (z)A 4 (z) + 2A-j (z)B^ (z) - A-j (z)A 3 (z)B 2 fz) + A^(z)B 1 (z) 

- A^z^ (z)] + N 2 (A>n,c|))[B 2 (z} - A 2 (z)B 2 (z) - 2B^ (z)B 2 (z)A 3 (z) 

+ A 2 (z)B 2 (z) + 2A 2 (z)A 4 (z)B 1 (z)B 2 (z) + 2A 1 (z)B 2 (z) 

- 2A 1 (z)A 3 (z)B 1 (z)B 2 (z) + A 2 (z)B 2 (z) - A 2 (z)B 2 (z)] 


A-j (z) 
A 2 (z) 


A 3 (z) 

a 4 ( z ) 

B-,(z) 

B 2 (z) 



The components of A(z) are determined as follows: 



[J G $4 + 2K p J G s3 + (K p + 2 K i J G ) $ 2 + 2K p K I s + KjKJy( s2 + K 3 s 
+ K 2 ) 2 - K 2 ] 


(2-17) 

(2-18) 

(2-19) 

( 2 - 20 ) 

(2-21) 

( 2 - 22 ) 

(2-23) 

(2-24) 
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The expressions for the discrete describing function of the CMG 
nonlinearity, N(A,n,(f>), may be found in [2]. In the present case, 
since the two axes are considered to be identical, and the couplings 
are symmetric, it is assumed that the amplitudes of the input signals 
at the two nonlinearities are identical. 

The stability equation is obtained by setting A(z) of Eq. (2-17) 
to zero. It is apparent that A(z) - 0 is of the form: 

A(z) = 1 + G ft (z) + Gg(z )N(A,n,<f>) + G c (z)N 2 (A,n,(t>) = 0 (2-25) 


which can be written as 

1 + Gg(T,n)N(A*n,<|>) + G c (T,n)N 2 (A,n,<i>) - 0 (2-26) 

by dividing both sides of the equation by 1 + G^(z), and setting 


z 


e j27r/n 


(2-27) 


Equation (2-26) shows that in the sampled-data system, there are 
four variable parameters in the stability equation in T, n. A, and <J>, 
whereas in the continuous-data system there are only two variables in 
A and w. 


2-3 . Exact Solution of the Stability Equation of the Two-Axis Sampled- 
Data LST System by Numerical -Iterative Techniques 


The numerical-iterative method described in Chapter 1 is now used 
to solve the stability equation of the two-axis sampled-data LST system, 
Eq. (2-26). 

In Eq. (2-26) we define 
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GgtT ,n) G^i + jGj-j 

(2-28) 

G c (T,n) = G r2 + jG I2 

(2-29) 

N(A,n,<f>) = N r1 + jNji 

(2-30) 

N 2 {A,n,(j>) = N r2 + j N j ^ 

(2-31) 


where G r -j , G^, G R2 , G^, N r -j , , N R2 , and Nj 2 are all real quantities. 

When Eqs. (2-28) through (2-31) are substituted into Eq. (2-26), 
and after grouping the real and imaginary parts, we have, 

A - A r + jAj = 0 (2-32) 

where 

A R = 1 + G R1 N R1 " G I1 N I1 + G R2 N R2 ' G I2 N I2 = 0 

A I = G I1 N R1 + G R1 N 11 + G I2 N R2 + G R2 N I2 = 0 

These two expressions represent two equations in four variables 
in A, T, <f>, and n. It is necessary to assign values to two of these 
variables and solve for the remaining two using the numerical -iterative 
method. In general, it is practical to assign values to T and n, and 
Eqs. (2-33) and (2-34) are solved to yield A and <j>. Alternatively, we 
may assign values to n and 4>, and. solve for A and T. The latter case 
has been carried out in the current study. 

As in Chapter 1, we define 

T 

A 




(2-33) . 
(2-34) 


x 


(2-35) 
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(2-36) 


Then, Eqs. (2-33) and (2-34) can be written as 

F(x) = 0 (2-37) 

The Newton-type quadratically convergent numerical method described 
in Chapter 1 can now be directly applied to this two-variable system. 

The following system parameters are used: 

T^fq = 0.1 ft- lb 

H = 600 ft-lb-sec 

Jg = 2.1 ft-lb-sec 2 

K = 5758.35 
o 

K-j - 1371.02 
K = 216 ft- 1 b/rad/sec 

r 

Kj = 9700 ft-lb/rad 
J v = 10 5 ft-lb-sec 2 
Y = 1.38 x 10 7 

Three sets of coefficients of coupling are used for this study: 

(1) K 2 = 0.1 K, =0.1 

(2) K 2 = 1.0 K 3 = 0.1 

(3) K 2 = 5.0 K, = 0.1 



For each integral n, and a fixed value of <j), a solution of Eq. (2-37) 
is attempted with an initial guess on T and A. A convergent solution for 
T and A together with the assigned n and (J> describes a sustained oscillation 
in the system. One of the characteristics of the numerical method is that 
the solution depends on the initial conditions set by the analyst. As 
shown in the earlier studies, it is quite possible that a system may have 
a multiple number of solutions for the same set of n and 4 >. An exhaustive 
search can only be conducted by a systematic scan of all the possible 
combinations of the initial values in the iterative process, unless one 
has general information at the outset on the number and the proximity of 
the solutions. 

Figure 2-4 shows the plots of the two solutions of T versus A for 
K 2 = K 3 = 0.1, and for relative large values of n. The results are shown 
only for <j> * 0°. When all the values of <J> are used, the solutions of T 

and A trace out a loop in the T versus A domain. For large values of n, 

these loops are very small, as the solutions are relatively insensitive 
to n. As n decreases, these loops become more pronounced, as indicated 
in Fig. 2-5. As n decreases, one of the solutions disappears, as the 
amplitude of A becomes exceedingly small. 

Figure 2-6 gives the computer results for = K 3 = 0.1 'that 
correspond to one set of solutions for <j> = 0°. Figure 2-7 tabulates the 
other set of results for K 2 = - 0.1 with $ - 0°. Figure 2-8 gives the 

tabulations of results for K 2 = K 3 = 0 . 1 , n = 4, and variable 

Figures 2-9 and 2-10 give the results with K 2 = 1.0 and = 0.1. 

Figure 2-11 shows the solutions for K 2 = 5,0 and = 0.1 with 0 - 0°. 

In this case, the amplitude of oscillation decreases as n decreases, so 
that there is less chance for the low-frequency oscillations to occur. 
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The results indicate that harder coupling (larger values of K 0 and K^) 
tends to reduce the chance of self-sustained oscillations. 



Amp! itude ( xIO D ) 

Figure 2-4. Amplitude and frequency ( 2-rr/nT ) of self-sustained oscillations for various sampling 
periods in the discrete-data two-axis LSI system; y - 1.38x10^. 





(sec) 
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LST SYSTEM-NUMERICAL SDLUTIDN OF DISCRETE DESCRIBING -FUNCTION • 
TWO Ha I S-EQUhL AMFL I TUBE CASE, K£= 1. GOOD- 01, K3= 1.0 GOD- 01 
GAMMA* 1 .3:30 OOP 0? H= 100 PHI= 0.0 


ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0.3 3 £ 0 0 0 0 0 0 0 D - 0 I 

0.6000 0 0 0 0 0 0 D - 0 6 

1 

0 . 334 £'££9 £5 0D- 0 1 

0.59 5 06793 0 70- 06 

o 

0.33 5 9 7 3 6 0 £ 3 D - 0 1 

0.5 i 4 ( 9 r 644 9 0 — ut> 

3 

0. 3366351 5 35 D- 01 

0 . 56685 1071 30- 06 

4 

0 . 33695665 120-0 1 

0 , 56£53£9 04 30-06 

Cj 

0 . 33695378 1 5D- 0 1 

0 . 56874 0569 1 D- 06 

6 

0 . 3 3 69 5986340- 0 1 

0 .54350 3 0 £ 9 0 0 - 0 6 

? 

0 . 3364 0 4 66 £ £ D - 0 1 

0 . 557 1 6994 39B- 06 

8 

0 .33676:37 1 £3P- 0 1 

0 . 5598 1 £54540- 06 

9 

0 . 33 6 985539 £ D — 0 1 

0 . 56 09o G5oo6B— 06 

10 

0.33 66. 94 06 0 4 D — 0 1 

0 . 56394 £754 UD— 06 

11 

0 . 33698 03£ £ 3 D - 0 1 

0 . 563 09639 1 SD- 06 


NB= 1 SAMPLING PERIOD* 3 .36972D-0SSEC ■ AMPLITUDE* 5.63342D-07 NT* 1 1 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE, K£= 1.0000-01, K3= 1.0 GOD- 01 
GAMMA* 1.33 GOOD 07 N= 36 PHI* 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 


0 0 . 336373 034 3D- 0 1 

1 0. 3474 33 153 3 D- 01 

£ 0 . 35 0£ 3754 0 £ D - 0 1 

3 0.351 1518195D-G1 

4 0 - 35 1574 33 6 3 D - 0 1 

5 0.351 3 OS 1 455D- 0 1 

6 0.351536746 6 D- 0 1 

7 0.351 S 6 354 6 £B - 0 1 


0 . 5 6 33 4 £3 5 5 0 D - 06 
0.544 3 3 533 3 3 D - 0 6 
0.5470 33 7 ££ 1 B - 06 
0.55123 0 31 1 7D - 06 
0 . 5 5 3 6 0 3 3 7 3 3 D - 06 
0 .554634361 70- 06 
0 . 5 56893 1 373D- 06 
U ■ 55665 05 4 4 3D — 0^* 


NB= 1 SAMPLING PERIOD* 3 . 5 1 849D- 02SEC AMPLITUDE* 5.56673D-07 NT* 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE > K'£* 1.0G0D-01, K3= 1.0000-01 


GAMMA* 1 ,330000 07 N= 92 
ITERATION SAMPLING PERIOD 

0 0 . 35 1 3498 1 ££D- 0 1 

1 0 . 3 6 3 £ 6 3 3 7 3 9 D - 0 1 

£ 0.3 6 6 3 3 £ 3 0 1 6 D - 0 1 

3 0 . 3673 53 4 75 9 D - 0 ] 


PHI* 0.0 

AMPLITUDE 
0 . 5 5 6679 033 1 D~ 0 6 
0.5 ^»Gt'4 cl' 4 1- r D — i.!6, 
0 . 5 39 j 7 055 7 3 D - ij 6 
0. 54 3554 3563 D- 06 
5D-06 


1 3£ 0D- 06 


4 

0.3678894 1 9 OB- 0 1 

0 .5458:59 

5 

0 . 36 8 0 89 74 0 7 D - 0 1 

0.54694 1 

6. 

0 . 3675 9 0 04 0 £B - 0 1 

0 .5494 03 

"? 

0 . 3679 86 6 0 4 4 D - G 1 

0 .549885 

3 

0 . 36.8 019178 6. B — 0 1 

0 . 56 1 34 1 

9 

U « •.> 6 o 1 y *7' *f' D " U 1 

0 . 55 06 1 ? 

1 0 

0. 36845 13633D- oi 

0.549525 

1 1 

0.368 3 £838860- 0 1 

0.549131 

18 

0 . 6 9 8 4 o 6 9 9 3 D - 0 1 

0 . 534625 

13 

0 . 369 t 5434 53D- 0 1 

0 ,53789 ,: i 

14 

0 . 3684 6533.3 00- 0 1 

0 .543968 

15 

0 .36?£7£5£04D- 01 

0 . 55 1 £ 04 


£539 D— 06 
5696 D~ 06 
4 63 3D- 06 
££4£D- 06 
8800D- 06 
7322D-06 


HD* 1 


SAMPLING PERIOD* 3 .67£73D-0£SEC AMPLITUDE* 


1 £050-07 NT* 15 


Figure 2-6a 
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LST SYSTEM-NUMERICAL SOLUTION 
Tim AX IS -EQUAL amplitude case 
l3BMMft= 1 .38 0 0 OD 0? N= 88 
ITERATION SAMPLING PERIOD 

0 0 .3678785804B-G1 

1 0.3804*510? ID- 01 

£ 0 . 384 05 78 5 74 D - 0 1 

3 0 . 3y 5 S 4 69 4 4 4 D — 01 

4 0.33 58 0 0(i 819 D- 01 

5 0 . 386 09 78 8 98 D ~ 0 1 

8 0 ■ 6 0 9 U ( ij 1 

7 0 . 3866869 363 B- 0 1 

8 0 . 3388538 39 OB- 0 1 

3 0 . 388 49485 84 0-01 

1 0 0. 366544 3397D- 01 

11 0 . 38658876830- 0 1 


DF DISCRETE DESCRIBING FUNCTION 
K 3= 1 .00 0 K -01 ’ K 3= 1.00 0 0-01 
PH I = 0.0 

AMPLITUDE 
0 .5518 047388 D- 06 
0 . 5876558668D- 06 
0 . 5897954 747D- 06 
0 .5 3394 05 78 3D- 08 
0 .5388538 0 5 9 D - 0 8 
Ci. 537348 1539D- 08 
0 . 5385987 083D- 08 
0 . 5 3 7871 8 8 8 4 D - 0 8 
0 . 5 4 0 7750 86 7 D - 0 6 
0 .537688678811- 08 
0 . 5 38997847 0 D- 08 
0 . 5 387388 06 90- 08 


NB= 


SAMPLING PERIOD* 3 .S6637D-02SEC AMPLITUDE* 5 .38489D-07 NT=11 


LST 

SY 

STEM- 

NUMERICAL SOLUTION 

DF DISCRETE DESCRIBING FUNCTION 

TUG 

AX 

I S' -EQUAL AMPLITUDE CASE 

. K£= i . 0 0 OB- 0 1 ■ K3= 1 .0000-01 

GAMMA* 

1 .38 

GOOD 07 N= 84 

PH I = 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 



0 

(i . 38 66 36- 64c ID — i.i 1 

0 . 5 3 6 4 8 8 5 3 8 1 D — 0 6 



1 

0.40 08 3 6 4 9 7 6 D - 0 1 

0 . 5 1 5 0 84 9 7 4 5 D- 06 



p 

0 .40391 00 1340- ul 

0.51 8 07 07747D- 06 



3 

0.405628 3 92 0 D - 0 1 ‘ 

0 . 582 1 6 092 5 ?D - 06 



4 

0 . 4 05S4 9 0375 B - 0 1 

0.52437174830-06 



5 

0.40618604450-01 

0 .5854 0833 1 9D- 06 



f. 

0 . 4 07 1 08423 1 0- 0 1 

0 . 5 8 3 9 8 74 4 3 9 1» - 0 6 




0 . 4 059663 U23D- 0 1 

0 . 5 3 0 9 9 984 7 6 D - 0 6 



o 

0 . 4 06373386 i 0- 0 1 

0 . 58643 0396 1 D- 06 



9 

0.4 06373654 8D- 01 

0 .58648 1 43750- 06 


ND= 1 SAMPLING PERIOD* 4 . 08373D-08SEC AMPLITUDE* 5.884140-07 NT* 9 


LST SYSTEM-NUMERICAL SOLUTION GF DISCRETE DESCRIBING FUNCTION 


TUG AXIS- 

EQUAL AMPLITUDE CAS 

GAMMA* 1 . 

3800 01* 07 H= 8 0 

I TERAT I DN SAM PL I NG PEP I GD 

0 

0.406373141811-01 

1 

0 .4819 079848D- 0 ! 

d 

0 .4861 34 95 1 8D- 0 1 

3 

0.4 8 7 66 0 75 0 6 1' - 0 1 

4 

0 . 4283895 084 D- 0 1 

tr 

0.42677904 1 ID- 01 

6 

0 .4294 1 6S292B- 0 1 

7 

0 . 469573957F.D- n l 

8 

0 . 48995995 1 0D- 0 1 

9 

0 . 42988 3384 ID- U 1 

10 

0.4 8 8 9 9745 2 9 D - 0 1 

11 

0.42972 34 9 2 0 D - 0 1 

18 

0.42123 585 8 8 B - 0 1 

1 3 

0.42 59 0 4 6 7 6 7 Ij - 0 1 

14 

0.4277 8 5 3 9 1 6 D - 0 1 

15 

0.48841 1 38 28 D- 01 


> K£= 1 . 00 0D- 01 > K3= 1.00 0D- 01 

PH I = 0.0 

AM PL I TUDE 
Ci . 5264 1 4 89 U41i— 06* 

0 . 5 0 1 4 7 0 4 4 78 D - 06 
0.5035481 161 D- 06 
Ci . 5 0 731 0 1 9 8 3 1' - 0 6 
0.5 0 9 3 6 8988 7 D - 0 6 
0 . 5 1 0389 1 74 3D- 06 
0 . 5 0974 1 1 8S7B- 0'6 
0 . 5 08488766SB- 06 
0 . 5 0 6 1356 3 7 7 D - 0 6 
0 .51806561 170-06 
0 . 5* 1 3 1564 0940- 06 
ft . 5 063 0364 9 00- 06 
0 .531 05038690-06 
0.51 86 98 6 4 1 8 0 - 06 
0.5147065691 B-06 
0.51374149590-06 


ND= 1 SAMPLING PERIOD* 4 .864 18D-08SEC AMPLITUDE* 5.137410-07 NT* 15 


Figure 2-6b 
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LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .380000 ft? M~ 7 A 
I TERAT I ON 5 HMPL I NG PER I DD 

0 0 .48841 1S328D- 01 

1 0.4464 US-1 67D -01 

2 0 .45 134 04 6*3 0D- 01 

3 0.46314 5 6 9 S3 D- 0 J 

4 0.464 00972690-Gl 

5 0 . 464463024 Op- 0 1 

6 0.466015 85 750-01 

7 0.455 0 3 62 1 6 71' - 0 1 


OF DISCRETE DESCRIBING FUNCTION 
K'£= l . 0 0 OP- 0 1 » K3- 1 . 0 0 0D- 0 1 
PH I = 0.0 

AM PL I TUBE 
0.5 1374 143590-06 
0 .4 £4 1 *7*4 c" 6 6 9 0— 06* 

0 . 4 3 5 0 S 0 u 6 2 3 0 - 0 6 
0 .4333325 0270- 06 
0.4902121 1730-06 
0.491 056 l 5020-06 
0.49 03 03 1 32 00-06 
0,49 06 1 0363 OP- 06 


NO- 1 SAMPLING PEP I OP= 4 .66 0560-02SEC AMPLITUDE* 4.9 05 030-07 NT- ? 


LST SYSTEM-NUMERICAL SOLUTION GF DISCRETE DESCRIBING FUNCTION 


TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .330000 07 N= 72 
ITERATION SAMPLING PERIOD 

0 0.455 0 66 4 3 7 3 P - 0 1 

1 0. 474349263c D-Ul 

2 0.4 3 03 2 7 0 73 6 It - 0 1 

3 0.43246901670-01 

4 0.4336 1634 4 lit- 01 

5 0.434 1 J 7 1 0 0 2 3 It - 0 1 

6 0.4 24 5 S 36 02 70 - 0 1 

7 0.434 5 9 1 33 2 3 0 - 0 1 


K£= 1 . 0 0 00-0 1 .♦ K3= 1.00 OP- 01 
PHI- 0.0 

AMPLITUDE 
0 .49050772 020-06 
0 .46 0267 6 06 Lip— U6 
0.46 It 0523230- 06 
U . 46*yy , c , c '"‘4*4 4 0 — u6 
0 .4654279 1 020- 06 
0 • 4 6* 6 1 U 0 y 9 4 6 D — U 6* 

U .466 U 1 c‘6 2 4 2 P — 0 6 
0 .4659959617 0 - 06 


HP* 1 


SAMPLING PER 100= 4 .345940- 02 SEC AMPLITUDE* 4.659330-07 


NT* 7 


LST SYSTEM-NUMERICAL SOLUTION DF DISCRETE. DESCRIBING FUNCTION 


TWO AXIS- 

-EQUAL AMPLITUDE CASE 

> 1 2= 1 . 0 0 0 D -01. K 2 * l.O 0 n p - n l 

GAMMA* 1 

.380000 07 H= 63 

PHI* 0.0 

I TERAT I ON SAMPL I MG PER I OP 

AM PL I TUDE 

0 

0 .464594232811-01 

0.46 5 9 8 3 81940-0 6 

1 

0 . 5 077 1 288 1 9P~ 0 1 

0 ,43055412 2 0 D - 0 6. 

2 

0.51414962770-01 

0 . 4 3 0 0 3 1 6 4670-06 

3 

0 .51678971 85 D— 0 1 

. 0.4 3 1 96 1 27330- 06 

4 

0.513 0 6 5 9 3 '=i 5 o - n l 

0.432 9 6 £ 3 £ 3 3 D - 0 6 

5 

0.518751991 00-01 

0 .43 3 3 7 3 3 3 0 3 0 - 0 6 

6 

, 0.519279 233 8 0 - 0 1 

0 . 43335556930— 06* . 

7 

0.51 9 8794 y 770—01 

0.43 3 3 3 3 1 2 8 5 D - 0 6 


ND= 1 SAMPLING PERIOD* 6 . 1 9280P-02SEC AMPLITUDE* 4 .233230- 07 


LST SYS TEN- HUME F I r A L SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .38 00 Op 07 N= 64 


ITERATION SAMPLING PERIOD 

0 0 .51927971420-01 

1 0 . 54677244950— 0 1 

2 0 . 55 44 0716 03 p - 0 1 

3 0.55767301 120-01 

4 0.55 9 3 2 3 7 2 1 -3 0-01 

5 0 . 56 0 1 946 04 0D- 0 1 

6 0 . 56 07441 7690- 0 1 


OF DISCRETE DESCRIBING FUNCTION 
j K£= 1 . 00 0D- 01 > K3= 1 . Ofi ft ji-n 1 
PHI* 0.0 

AMPLITUDE 

0 . 4 3 3 32 7 8 5 8 8 D - 0 6 . 

0 . 39 1 06 7289 OP- 06 
0 . 3 8 8 9 3 8 0 3 78 0 - 0 6 
0 .389650195311-06 
0 . 3 3 9 934 3 0 460- 0 6 
0 .83996*842 1 70“ U6 
0 , 3399899 1 2 6P- 06 


HO- 1 SAMPLING PERIOD* 5 .60744 D- 02SEC AMPLITUDE* 3.900100-07 


NT* 6 


Figure 2-6c 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE • K2= 1. 0000-01* K3= 1. 0000-01 
GhMMh= I . So 0 COD 0? N= 60 FHI = 0.0 

1 7 ERR T ION SAMPLING PERIOD AMPLITUDE 

0 0.56074407840-01 0.39001043010-06 

1 0.59438551840-0 ! 0 . 3377863 1 840- 06 

8 0.60 34 3573990- 0 1 0 .334 09655050-06 

3 0 .6 0756496 1 3D — 0 1 0 . 3 336976694 D— 06 

4 0.609689: 0 3 9 9 D - 0 1 0.33664517440-06 

5 0.61 03 1 486 06 D- 0 1 0 . 338 1 7 7434 7D- 06 

6 0.61 13998606D-01 0.3381931 0130-06 

7 0.61139934770-01 0.36890348500-06 

8 0.61091761 63 D - 0 1 0 . 3 36 1 3 56 964 0- 0 6 

9 0.61 136 16894D-01 0 .336446580-50-06 . 

10 0.611 36 064650- 0 1 0 . 33778 1 37860- 06 

1 1 0.61601868650-01 0 .33135846940-06 

12 0.6160 8975 190- 0 1 0 » 3 3 0 1 6 8 6 886 D- 06 

13 0.6161 3374360- 0 1 0 . 3£3?6?7 1 550- 06 

14 0,616181161 90- 01 0 . 369592 1 £640- 06 

ND= 1 SAMPLING PER IDD= 6 . 16648D-06SEC AMFL I TUOE= 3.693510-07 HT=14 


LST SYSTEM-NUMER I CRL SDLU IT ON OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE t K8= l.uOOD-Ol) K3= 1.0000-01 
GRMMA= 1.330000 07 M- 56 PHI = 0.0 

ITERRT ION SAMPLING PERIOD AMPLITUDE 

0 0 .6 1864 794690- 01 0 .36935098780- 06 

1 0 .654 33 01761 0 - 0 1 0 .66395775 03 D - 0 6 

£ 0.66451? 3 6460 — 01 0 . £6 08653 3860 — 06 

3 0.66 95 7 979 4 00 - 0 1 0 . £ 58 9 4 94 £ 4 4 0 - 06. 

4 0 .67866044 060- 0 1 0 . 65743746 050- 06 

5 0 .6737 01 98430-01 0 .25643975950- 06 

6 0.67438969030-01 0 .65688798060- 06 

7 0 . 674393884 70- 0 1 0 . £557 i 570660- 06 

8 0.6744092 38 8 0-01 0 .8559 0 9 4 3 3 3 D - 0 6 

ND= 1 SAMPLING PERIOD- 6 .744 17D-06SEC HMF‘LITUDE = £.553430-07 NT= 8 


LST SYSTEM-NUMERICAL SOLUTION OP DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* K6= 1.0000-01, \C3~ l.n 000-01 
GAMMA= 1.380000 07 N= 56 FHI= 0.0 

ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 .6744 1 669770- 0 1 0 .85584305790- 06 

1 0 .73635 077670- 0 1 0.1 596 1 66 0 050-06 

£ 0 .75 1 38699330-0 1 0.11 3 03763560- 06 

3 0.75441970060-01 0 . 1 4 1 4355664 0- 06 

4 0 . 7548893 1 3 70- 0 1 0.15 39 3 8 1 i 1 1 D - 0 6 

5 0.75400376340-01 0 . 15779657470- 06 

6 0 .7534 0563630- 01 0 . 16009714630-06 

7 0.75346758540- 0 1 0 . 1 5 8 9 l 1 3 3 34 0 - 0 6 

8 0.7656 36 1 7 1 3 0 -01 0.1 0 56 £ 6 8: 6 390-0 6 

9 0 . 76 096346760- 01 0.1 £9753733 1 0- 06 

10 0 .75756465640- 0 l 0 . 1 448 1 236980- 06 

11 0 . 755 1 35964 10-01 0.1 53 1 4 975930- 06 

16 0.75416131680-01 0 . 1 5787 0533 70- 06 

13 _ 0 .7535353354 0-01 0 . 1597 0089930-06 

14 0.75373430930-01 0.15790774930-06 

15 0.7 53:3 086 1 350- 0 1 0.1 6. 0 0 0554 090- 06 

ND= 1 SAMPLING PERIOD- 7 .533080- 08SEC HMPL I TUDE= 1.600060-0? NT=15 


Figure 2-6d 



LST SYSTEM-NUMERICAL SOLUTION 
TWO FIX IS -EQUAL HURL I T UDE CASE 
GRAMA- 1 .SCO 0 0 D 0 7 A = 1 0 0 
I TERRI I DM SAMF'L IMG. PERIOD 

0 0.15300 0 0 0 0 0 D - 0 1 

1 0. 15131 29565D- 01 

2 0.151 2363339D- 0 1 

3 0.150953849SD-01 

4 0.1507480998D-01 

5 0.1 5075774 19D- 01 

6 0.1503043353D-01 

7 0.150593019711-01 

3 0. 15059?S'776D-01 


OF DISCRETE DESCRIBING FUNCTION 
, K£= K3=0. 1 
PH I = 0 . iJ 

AMPLITUDE 
0.500 0 0 0 0 0 0 0 D - 0 5 
0 .512319453311-05 
0.51 33983252 P - 0 5 
0 . 52 1 88 08809 D-05 
0 .524 0422 024 D- 05 
0 . 523 1 344 0480- 05 
0.522576483 OP- 05 
0.52531 £32720- 05 
0 . 525356 225 7D- 05 


ND= 1 SAMPLING PERIOD- 1 .50598P-02SEC AMPLITUDE* 5. £53870-06 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .38 00 0 Ei 07 H» 96 
ITERATION SAMPLING PERIOD 

0 0. 1 50597? 177D- 01 

1 0 . 15246.8 64 7 3D - 01 

£ 0. 1545679 048P- 01 

3 0-1 556856 3 0 6 D — 0 1 

4 0 . 1562697461 D- 01 

5 0.1 565694849B- 0 1 

6 0 . 15661 06863D- 01 


OF DISCRETE DESCRIBING FUNCTION 
, K2= K3=0. 1 
PH I = 0.0 

AM PL I TUBE 
0 .5 £538654 9 8 D- 05 
0 .542421 0943D- 05 
0 . 535 068779 1 T U5 
0.53065470 9 9 D- 0 5 
0 . 528255 0 09 1 D- 05 
0.527 09 12439D-05 
0 .52686 1 2 07 4 D- 05 


ND= 1 SAMPLING PERIODs= 1 .56629D-02SEC AMPLITUDE* 5.26612D-06 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA= 1.3S0 0 0D 07 N= 92- 
ITERATION SAMPLING PERIOD 

0 0.1 566c*y9b'5SD- u 1 

1 0.1 5869283 07D- 01 

2 0.161 02695 10D- 01 

3 0.1 62268667 5 D - 0 1 

4 0 , 1629157791 D-01 

5 0.16325 0 032 6 D-01 

6 0 . 163308464 1 D-01 


OF DISCRETE DESCRIBING FUNCTION 

» K2= 0=0.1 
PH I = 0.0 

AMPLITUDE 
0.52661 16798D-05 
0 . 5447574 799D- 05 
0 . 53679983 1 4 D- 05 
0 . 5,32 04 0788 1 1'- 05 
0 .5294671 1 35D-05 
0 . 52 8 1 9 05 68 0 D- 0 5 
0 .527968 1922 D- 05 


ND= 1 SAMPLING FERIOD= 1 .63327D- 02SEC AMPLITUDE* 5 .27664 P- 06 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .38 GOOD 07 N= 88 
ITERATION SAMPLING PERIOD 

0 0.1 633272606D- 01 

1 0.1 655487224 D- 0 1 

£ 0.1 68 1 022992D- 0 1 

3 0 . 1694596 034 D- 0 1 

4 0.1701 6 55283D- 0 1 

5 0.17 053 1 0 04 0D- 0 1 

6 0 . 1 7 06 0596 1 6L- 0 1 


OF DISCRETE DESCRIBING FUNCTION 
* K£= K3=0. 1 
PH I = 0.0 

AMPLITUDE 
0 .527 6 6 36 6 8 3 D - 0 5 
0 .5469 1 592440-05 
0 .538474393 IB- 05 
0.533445651 1 D-05 
0 .5307374841 D- 05 
0 .5293798727 D-05 
0 .5291479667D-05 


ND= 1 SAMPLING PERIOD* 1 .70625D-02SEC AMPLITUDE* 5.28804D-06 


Figure 2-7a. 



LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING 
TWO AXIS-EQUAL AMPLITUDE CASE > K2=K3=0. 1 
GAMMA* 1.33 GOOD 07 N= 34 
ITERATION SAMPLING PERIOD 


FUNCTION 


0.170634 9 398 P - 0 1 
0 . 173031691 OP- 01 
0.1 75826433 OD- 01 
0.1 773 1 59629P- 0 1 
0. 1 78089 145 ID- 01 
0 . 173489897 OP- 01 
0. 1 7853572 1 SB- 01 


PH I = 0.0 

AMPLITUDE 
0 . 5288 0353 05P- 05 
0 .54 9 3 965 1 86 P — 05 
0 .54 03218271 P-05 
0 .534996 0G75D- 05 
0 .533 1 337 53 5 P- 05 
0 .530 6 9635220- 05 
0.53 043933 1 9P- 05 


NP= 1 


SAMPLING PERIOP= 1 . 7S605P-03SEC AMPLITUDE- 5.30055P Ob 


LST SYSTEM-NUMERICAL SOLUTION OP DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE » K2=K3=0.1 
GAMMR= 1 .38 0 OOP 07 N= 8 0 PHI= 0.0 


ITERATION SAMPLING PERIOD 

0 0 . 1 7 86 04 584 4 D - 0 1 

1 0.181 1982337D- 01 

2 0.1 8 4 8 9 2 2 8 y 9 D — 0 1 

3 0.1 85934 1 236 P- 0 1 

4 0 . 1 86 7 84 7 392 D ~ 0 1 

5 0.1 3722569S3P- 0 i 

6 0 . 1 873473958D- 0 1 


AMPLITUDE 
0.53 0 05534 1 6D- 05 
0 .5519351613D-05 
0 . 54237 09952D- 05 
0 . 5367 1 654640- 05 
0 . 53 3 69 64 7 7 9D- 05 
n . 532 161391 5D - 05 
0 .531862561 3D- 05 


NB= 1 SAMPLING PERIOD* 1 .87366D-02SEC AMPLITUDE- 5.31 44 OD 06 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AX I S-EQUAL AMF'L I TUDE CASE j K2= K3=0. 1 

PH I = 0.0 


GAMMA= 1 .38 00 OD 07 N= 76 
ITERATION SAMPLING PERIOD 

0 0.1 37365966 6D- u 1 

1 0 . 1 90 1 82 081 4 D- 0 1 

2 0.19361 16657D- 01 

3 0.1 9543 0 037SD~ 0 1 

4 0.1 9637 0387 OD- 0 1 

5 0. 196857552 ID- 01 

6 0 . 1 97 0 1 06 583 D- 0 1 


AMPLITUDE 
0 .5314404 6 8 5 D - 0 5 
0 .55487497270-05 
0 . 5446568 1 1 6D- 05 
0 . 5 3 8 6 3 7 5 0 3 £ D - 0 5 
n . 5354355 062D- 05 
0.53 38 01622 3 D-05 
0 .53344 121 32D-05 


NDa 1 SAMPLING PERIDD= 1 .97029D-03SEC AMPLITUDE* 


>29 8 7 D- 06 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE * K2=K3=0.1 
CAMMA= 1 .380 OOP 07 N= 73 PHI = 0.0 


ITERATION SAMPLING PERIOD 
n . 1 97 02898 1 7D- 0 1 
0 .20009833750-01 
0 .2 039 1 98 0821*- 0 1 
0.205944279*0-0 1 
0.206 93 935 9 OD - 0 1 
0 . 207530 08 1 9 It- 0 1 
0 . 20772 05288D- 01 


AMPLITUDE 
0 . 5329874 32 1 D- 05 
0 . 558 1 692582D- 05 
0 .547222869 1 D-05 
0.54 0 7? 6888 0 D-05 
0 . 5373 92 0 4 0 0 D - 05 
0 . 5 3 5 6 5 0 7 8 4 9 D - 0 5 
0 . 5352 094 1 85P- 05 


NP= 1 SAMPLING PER 100= 2 . 07739D-02SEC AMPLITUDE* 5.34734B-06 


Figure 2-7b. 
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LST SYSTEM-NUMEP I CAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA® 1 .38 00 OB 07 N= 68 
ITERATION SAMPLING PERIOD 

0 0.2077389646P-01 

1 0 .8 I 1 09y0716D~ 01 

2 0.21 538 1 9 1 75D- 0 1 

3 0.21764885040-01 

4 0.21581714400-01 

5 0.2 1842 05 08 ID- 01 

6 0 .2 196545462D-01 

7 0.218673871 ID- 01 


OF DISCFETE DESCRIBING FUNCTION 
K£= K3=0. 1 
PH I = 0.0 

RMPL ITUDE 
0 . 534 733829 0D- 05 
0 .56 18S43722D-05 
0 . jjOI 24359 7 0- 05 
0.54 324 32679 D- 05 
0 .53961 12375P-G5 
0 .5377523 077 D- 05 
0. 5372 1 2 051 9 p- 05 
0 .5367283776D- 05 


NB= 1 SAMPLING PERIOD® 2 . 1 9698P-02SEC AMPLITUDE® 5.36633D-06 h- 7 


LST 3YSTEM-NUMEPICRL SOLUTION 
TWO AXIS-EQUAL RMPL I TUBE CRSE 
GAMMA® 1 . 38 0 0 0 D 07 N= 64 
ITERRTION SAMPLING PERIOD 

0 0.219697 9 44 4 P - 0 1 

1 0 . 2233797565D- 0 1 

2 0.2 £32 06 5 S 1 7 P— 0 1 

3 0.230759 001 5D- 01 

4 0 . 232 07219 0 0 D- 0 1 

5 0 . 2 327 4 89 02 9 D- 0 1 

6* 0.2330 3 2579 3 D - 0 1 

7 0.233 054 121 6 D- 01 


OF DISCRETE DESCRIBING FUNCTION 
, K£=K3=0.1 
PHI= 0.0 

RMPL I TUBE 
0 . 5366325 1 66Ii- 05 
ij . 566 U864357B— 05 
0 . 5534 1 696830- 05 
0 .54603235160-05 
0. 54214830970-05 
0 . 54 0 1 6 1 8738D- 05 
0 . 5 39 5 03732 2 D - 0 5 
0 . 539 033 1 95 1 D- 05 


HD® 1 SAMPLING PERIOD® £ . 33080D-02SEC RMPL I TUBE = 5 .38 932 D- 06 N= 


LST S Y STEM- N UME R I C fiL SOLUTION 
TWO AX IS- EQUAL RMPL ITUDE CASE 
GAMMA® 1 .380 0 011 07 H= 6 0 
ITERATION SRMPLING PERIOD 

0 0 . £ 3 0 3 0 0 8 7 6 D - 0 1 

1 0 • £37 1 546467D— U 1 

2 0.242 6 3 8 68 6 9 P - 0 1 

3 0 . 24 5 53 573 89D - 0 1 

4 0.247023 63 92 D - 0 1 

5 0.247738755 ID- 01 

6 0.2481 500 177D-01 

7 0.248 156 038 ID- 01 


OF DISCRETE DESCRIBING FUNCTION 
. K2= K3=Q. 1 
PH 1= 0.0 

RMPL ITUDE 
0 . 53893 1 7 022D- 05 
0.5709270125D-05 
0 . 557287390 OP- 05 
0.5492 6 2472 3 D - 0 5 
0.54 5 0 6 94 726 D - 05 
0. 5429575 152D- 05 
0. 5421 7853 19D- 05 
0.541 7289555 D- 05 


ND= 1 SRMPLING PERIOD® 2 .431S5D-02SEC RMPL I TUBE® 5.4161 ID-06 N= 7 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO R: % r S- EQUAL AMPLITUDE CASE, K£®K3=0.1 


GAMMA® 1 .38 

0 0 OB- 07 

H= 56 

PHI= 0.0 

ITERATION 

SRMPLING PERIOD 

AMPLITUDE 

0 

0 . 848 18457751'- 0 1 

0 . 54 161 09 1 9 OB- 05 

1 

0 . 252722259 OB- 0 1 

0.57 6 5 3 4 1 9 8 3D- 0 5 

2 

0 , 259 0 0 371 1 6 P - 0 1 

0 . 56 1661 6*8 1 CD— U 5 

3 

0 .26231 8U446D- U 1 

Ci . 5 5 3 0 3 9 5 0 4 3 D - 0 s 

4 

0.264 01 7762 ID- 01 

. 0.548528 14 70D- 05 

5 

0.26 4 83 9 0472 D- 01 

0.54 6 24345 3 1 1' - 0 5 

6 

0 . 2652964524 D- 01 

0,54 5 32 8 6 1 4 6 D - 0 5 

7 

0 .26533 0" 

’74 ID- 01 

0 . 54492 1 1 264 D- 05 


HD® 1 SAMPLING PERIOD® £ .65363D-02SEC RMPLITUDE= 5.44773D-06 H® 7 


Figure 2 -7c. 



LST SYSTEM-NUMER I CFtL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA® 1.38 GOOD 07 M~ 52 
ITERATION SAMPLING PERIOD 

0 0 . £6536£6£94 D-0 1 

1 0.2704516365D-01 

£ 0.27771£68£9D-01 

3 0.281538013211-01 

4 0.2834965647D-01 

5 0. 284497587811- 01 

6 0.284981 7 152D- 01 

7 0. 2350309567H- 01 


OF DISCRETE DESCRIBING FUNCTION 
* K2» K3=0. 1 
PHI® 0.0 

AMPL I TUBE 
0 .54477260570-05 
0 .5331 0664740-05 
0.56 6 8 9 £ 9 3 0 5 It - 0 5 
0 .557520150 3 D - 0 5 
0. 552644 016SD- 05 
0.55 01 6533 7 ID- 05 
0 .5491 072646D-G5 
0.5437526844D-05 


NB® 1 SAMPLING PERIOD® 2 .S5067D-02SEC AMPLITUDE- 5.4S557D-06 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* K2®K3*0.1 


GAMMA® 1.3 3 

GOOD 07 N= 48 

PHI ® 0.0 

ITERATION 

SAMPLING PERIOD 

AM PL I TUBE 

0 

0 . 235 06653 8 2D- 0 1 

0 . 548556876SB- 05 

1 

0 . 29 0 8 1 32 1 9 0 D - 0 1 

0 . 59 o 91954 09 D - 05 

£ 

0 • £ d ’? y 7* & S 3 Ii — u 1 

0.5731 615644 D- 05 

3 

0.3 0375369 1 6D- 0 1 

0 . 5 G 2 9 2 7 8 3 ( > •' D — u 5 

4 

0.306 0 375 8 2 8 D - 0 1 

0 .557 62247 640— 05 

5 

0 . 3 071 9375 1 7D- 0 1 

0 . 554934S977B- 05 

r> 

0 . 3 077733 0 74 D- 01 

0.5 53 7263 22 CD -05 

7 

0.3 073439653 D- 01 

0 .553425253 1 B- 05 


NIi= 1 SAMPLING PERIOD® 3 . 07S88D-02SEC AMPLITUDE- 5.53163D-06 


TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1 .38 GOOD 07 N= 44 
ITERATION SAMPLING PERIOD 

0 0.3 0733771 4 ID- 01 

1 0 .3144439 053D- 01 

£ 0.32447 0936 ID- 01 

3 0 . 3397263 07 0D- 0 1 

4 0 - 3324 066935D- 0 1 

5 0.33576 S' 0 36 S D -01 

6 0.3344 5 1298 ID- 01 

7 0 . 3345 72 7 029 D - 0 1 


DISCRETE DESCRIBING FUNCTION 
K2= K3=0 - 1 
PH I = 0.0 

AMPLITUDE 
0 .5531631 338D- 05 
0 . 6 0 03643 0 1 GO— 05 
0 . 53 08 1 558 04 P- 05 
0.56959141 73D-05 
0.5637937 05211-05 
0 . 56. 0663 1 697 D - 05 
0.5595031 02 OB- 05 
0 ,55923618330-05 


LST SYSTEM-NUMERICAL SOLUTION OF 


ND= 1 SAMPLING PERIOD® 3 .34613D-02SEC AMPLITUDE® 5.53839D-06 


Figure 2-7d 
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LST SYSTEM-NUMERICAL SOLUTION 
Tun AXIS-EQUAL AMF L. I TUPE COSE 
GANNA = 1 . 38 OH OP 01 * N- 20 
I i L Rfl 1 i UN I -f'lMPl. 1 MG PER 1 Cl D 

0 0 . 6 0 0 0 0 0 0 0 0 0 D - 0 t 

1 0 . 6480379984 D- 0 1 

2 Q . 6 6 78 026 u i 9 D - 0 1 . 

3 0. 67695821 2 Ofl- 01 

4 0. 6S 13 04334 lit- 01 

5 0 .683425 0428 D- 0 1 

6 0.6345 08 3594 D- 0 1 

7 0 .68 44 64 554 3D- 0 1 

3 0 .6346078 1 7-JD-01 


OF DISCRETE DESCRIBING FUNCTION 
K2«K>0. 1 
PH I - d.U 

AMPL. I IUDE 
0.10 0 0 0 0 0 0 0 0 D- 04 
i.i ,761 1 6c'9y54D- 05 
0.71 72493 1 89 D- 05 
0.69968221 8 8 D- 05 
0 .69184331 5 0D- 05 
0 , 688253 1 1 39D- 05 
0 . 6863 1 9 06 1 2D- 05 
0 . 6 8 5 8 1 6 2 9 0 8 D - 0 5 
0.685474493 0D-05 


ND= 1 SAMPLING PERIDD= 6 .84679D-02SEC AMPLITUDE- 6.85169D-06 


N= 8 


LST SYSTEM- NUME PI CAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1.5 8 0 0 0 D 07 N = 84 

ITERATION SAMPLING PERIOD 

0 0 . 6 3 467379040- 01 

1 0.678 9 6 8 8 9 1 j 3 D - 0 1 
0.6- 38 7 38 1 1- 58 B — IJ 1 

3 0 . 583 1 5 92 1 93 D- 0 1 

4 0 .57251 47888 D-0 1 

5 0. 5 73 69725 79 B- 01 

6 0 . 57923S3S75P- 0 1 

7 0.5 8 2 0889924 D - 0 1 

8 0.5835654 75 ID- 01 

9 0.5 8 4262 9 6 5 2 D - 0 1 

10 0 . 584287 0822D- 0 1 


DF DISCRETE DESCRIBING FUNCTION 
K2=K3=0. 1 
PHI = 0.0 

AMPLITUDE 

0 . 6 8 5 1 6 9 2 6 6 3 0 - 0 5 
0 . 3849797353B-05 
0 . 4653644434 D— 05 
0.51 1355551 7B- 05 
0.61 365 08 1 63 D- 05 ‘ 

0 . 64878854 5 3D— 05 
0 . 6 4 3 5 5 S2 S 88 B - O 5 
0 .64 06344573D- 05 
0 . 6 3 9 0 8 Ci 3 tj 9 6 D - 0 5 
0 . 638568 1 095D- 05 
0 .638 1 238265D- 05 


HD- 1 SAMPLING PERIOD- 5 .84322D-02SEC AMPLITUDE* 6.381 06D- 06 N=10 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1 .33 00 OB 07 N= 28 


ITERATION 

SAMPLING PERIOD 

0 

0 . 58432 1 552 9 D - 0 1 

1 

0 . 532 0296 04 7D - 0 2 

C ' 

0 . 5 5 0 38 04 67 7 D - 0 1 

3 

0 .531 82423 04D- 0 1 

4 

0.521 021914311-01 

5 

0.51 5285290 3D- 0 1 

6 

0.51243443 0 8 D - 0 1 

7 

0.51 1 174 151 7D- 01 

c* 

0.51 06 04 6578 D— 0 1 

. 9 

0.50 9 5 3 0 2 3 7 9 D - 0 1 

10 

0. 5 09 56841 64 D- 01 

11 

0 . 5 095 12 04430-01 

12 

0.50 96 9 3 7 82 5 0 - 0 1 

13 

0 . 5 09694 0058D- 0 1 


OF DISCRETE DESCRIBING FUNCTION 
K2*K3=0. 1 
F'H I ~ 0.0 

AMPLITUDE 

0 . 638 1 063 i 2 1 D- 05 
0 . 4 07936 2 2 34 1> - 05 
0.4770 3 8 8 7 8 9 D - 0 5 
0.529893 735 ID- 05 
0 . 565 0 4 672 £ 7 D - 0 5 
0 .58589 1 29690-05 
ij .59726 04 748D— 05 
0 . 6 0 3 0 0 5 2 9 0 9 D - 0 5 
0 . 6 0 563 0 1 3 8 8 D - 0 5 
0.61 009371 01D-05 
0.60 3 8 8 7 3 9 4 2 D - 0 5 
0 .6074034692D- 05 
0 . 6 0 8 0 5 017 0 0 D - 0 5 
0 . 6 077955 1 58 D- 05 


ND= 1 SAMPLING PERIOD- 5 . 09695D- 02SEC AMPLITUDE- 6.07377D-06 H=13 


Figure 2-7e 



lst system-numerical solution df discrete describing function 

TWO AXIS-EQUAL AMPLITUDE CASE K2=K3=0.1 


GAMMA= 1 .33 

GOOD 07 N- 38 

PH I = 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0 - 5 0 96 9 4 53 0 1 D - 0 1 

0 . 6 0 7 S' 7 7 OS 6 S D - 0 5 

1 

0.51084 9 7 98 8 D - 0 1 

0.411916414 0D- 05 

£ 

0 .43434788840- 0 1 

0 . 47570 1901 6D- 05 

3 

0 .46396836 08D- 0 1 

0 . 58868346890— 05 

4 

0. 46 0457 38 S3 D~ 01 

0 . 5587945 1 860-05 

c 

0 . 4 56 08 0 8 1 05 D- 0 1 

. 0 . 5 7 088 1 0 94 7 D - 05 

6 

0.45333561 130- 01 

0.5795 1556580-05 

7 

0.45879818 9 9 D - 0 1 

0 .584 155 184 3D- 05 

8 

0.458194466 5 D ~ 0 1 

0 . 5 866 8 1 8 5 85 D - 05 

9 

0 .45473 008790- 0 1 

0 . 5 7 8’ 7 9 1 6 8 3 8 D - 0 5 

10 

0.4 5c63 1 56 UoD- 0 1 

0 . 53 1 4 044 1 94 D- 05 

11 

0 . 4 53 7488 954 D - 0 1 

0 . 5854 1 4 7556D- 05 

IE 

0,4584681 04 3 D - 0 1 

0 .588 1 1 98 35 3D- 05 

13 

0 .43644 131 64 D- 01 

0 . 56 087 0987 0D- 05 

14 

0.463 3 4 6536 3 D - 0 1 

0 . 58 87 431 6 1 6 D - 0 5 

15 

0.4573 083 7 7 0 D - 0 1 

0 .537 044:39 1 1 p- 05 

16 

0 .45457 3 S 3 8 9 L - 0 1 

0 . 58 7 3 7 3 8 7 5 7 D - 0 5 

17 

0 . 4 53 04 £30770- 0 1 

0 , 5 8 8 3 699 8 790- 0 5 

13 

0 .4519353759H-0 1 

0 . 5 9 0 8 39 07570-05 

19 

0 .451S646046D-01 

0 . 5 8 8 5 84506 9 O — 0 5 

£0 

0 .451875651 3P- 01 

0 . 5 8 788 8 0 0 9 6 O - 0 5 


ND= 1 SAMPLING F'ERIOD= *4 .51876D-08SEC AM PL I TUBE = 5 .878880- 06 

LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE K2=K3-0.1 


GAMMA = 1 . 38 

Liu OD 07 N= 36 

PHI - 0.0 

ITERATION 

sampling p epic o 

AMPLITUDE 

0 

0.451875651 30-01 

0 . 5 8 788 8 0 0 9 6 0 - 0 5 

1 

0.45434418050- u 1 

0 . 4 0 795397 34 0 - 0 5 

2 

0.4 3 194 98 1350- 01 

0 .47048 1864 7D- 05 

3 

0.41 875353370- 0 1 

0 .51 587535480- 05 

4 

0.41161 89664 0- 0 1 

0 .543 3 8 5 8 0 3 0 1' - U 5 


0 .4 0794 0 0964 D- 0 1 

0 . 5 59 348 3 3760- 0 5 

6 

0.406151 68S8D- 0 1 

0 . 56 78 09 1 4 3 00— 05 

r y 

0 . 4 054 5806470- 0 1 

0.571 9556 0 1 30- 05 


0.40517354010-01 

0.5 7384 374 570- 05 

*3 

0.40580845140-01 

0.57334671740-05 


HD= 1 SAMPLING PERIOD- 4 . 05884D-08SEC AMPLITUDE^ 5 .789360- 06 


Figure 2-7f. 
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LST SYSTEM- MUM EP I CFIL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EOUAL ' AMPLITUDE CASE K2=K3=0.1 


GAMMA = 1.3G 

0 0 OP Or N= 4 0 

PH I = 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0.40528417940-01 

0 .57293 59 1 39 D - 0 5 

1 

0.41 0537576 6D- 01 

0 . 3 9 9 3 9 4 5 0 9 9 D - 0 5 

2 

0 . 33864 07364 1>- 0 1 

0. 46425 147S5D- 05 

3 

0 . 375496435 OB- 0 1 

0 . 5 096797 06.9D- 05 

4 

0 . 379439922 OP- 0 1 

0 . 54853 04 1 OSD- 05 

5 

0 . 37266 029 1 1 V~ 0 1 

n , 55962465 07D- 05 

6 

0 ■ 369664 72c 9P- 0 1 

0 .562 9 6 38 S 7 6 D — 0 5 

7 

0 ■ 366’ 1.466579P- M 1 

0 .564 6 6 0 1 £ 6 7 D - 0 5 

3 

0.3671 742546D- 01 

0 .56690536U7D- 05 

9 

0. 367 I 836255 D- 01 

0 . 563 776382 1 0- 05 

1 0 

0.3674581 1 25P- 0 1 

0 .56 1 3664939 D- 05 

11 

0 . 36 712907 7.9 B - 0 1 

0 . 5 6 3 5 2 8 1 0 £ 9 D - 0 5 

12 

0. 3674 02271 8 D- 01 

0 .561 6 0 5231 8 D - 05 

13 

0 .367 U94 069 UO- U 1 

0 . 56366844690- 05 

14 

0.367241 041 3D- 01 

0 .5624874991 D- 05 

15 

0 . 3 669 1 7c62 1 D — 0 1 

0 . 56 4 5 1 32 7 7 8 D - 05 

16 

0 . 36 6 9 2576 6 4 1' - 0 1 

0.5 6 4 0 1 3 0 3 4 8 B - 0 5 


ND- 1 SAMPLING PERIOD* 3 .66940D-02SEC AMPLITUDE* 5.63S92P-06 N=16 


LST SYSTEM-NUMERICAL SOLUTION. OF DISCRETE DESCRIBING -FUNCTION 


TWO AXIS-EOUAL AMPLITUDE- CASE K2=K3=0.1 


HNflH- 1 .38 

0 0 UD 07 N= 44 

PHI = 0.0 

HERAT inn 

SAMPLING PER I CD 

AMPLITUDE 

0 

0 . 36694 04 0 1 9 p - ij l 

0 . 563 89 1997 6 D - 0 5 

1 

0 .374 0 1 3853 0 0 - 0 1 

0 . 39 089 1 55330- 05 

a 

0.349101053211-01 

0 .45 3 8 9 3 4 6 3 4 0 — 0 5 


0 .3431 1 36-753 D- 0 1 

0 .53367 0882 ID- 05 

4 

0.341 1 1959 ODD- 0 1 

0 . 54834986 04 D— 05 

cr 

0 .3379425697D- 0 t 

0 . 5 54 1 82 0 3 0 o 0 — 0 5 

£. 

0 . 33638:2938 ID - 0 1 

0 . 55682 1 1 9720- 05 

7 

0. 3354 3674 7 OB- 01 

0 . 559342' 1 7 040- 05 

8 

0 . 3354678742D- 0 1 

0 .556 0 5 3 3 5740- 0 5 

9 

0 . 336 987893 ID- 0 1 

0 . 54 57 1 1 02 1 40- 05 

10 

0 . 335804236 1 p— 0 l 

n . 558477554 1 p- 05 

1 1 

0 . 3354 1 43435D- 0 1 

0 .5555 8 6 3 0 0 9 0 - 0 5 

12 

0 . 33 4 8 7 6 6 ? 3 6 0 - 0 1 

0 .5539398 1 8 00— 05 

13 

0.334375355711-01 

0 . 56?!-:66:-; nS6D- n5 

14 

0 . 33527 1 9 08 1 P- 0 1 

0 .55263483490-05 

15 

0 .33523376860- 01 

0 . 5 5 9 696 6 8 6 4 0 — 0 5 

16 

0 .33521 1 9736 D- 01 

0 .557 5 0 7 0 0990- 0 5 

1? 

0.33522247 130-01 

0 . 556 6 0 197 5 8 P - 0 5 

18 

0 .33524879970-0 1 

0 . F, F, 6 3 4 8 £ 6 3 0 - 0 S 


HD= 1 SAMPLING PERIOD* 3 .35304B-02SEC AMPLITUDE* 5. 56 04 ID- 06 N-18 


Figure 2-7g 
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LN SI 


LSI SYSTEM -NUMERICAL SOLUTION 
TWO AXIS-EQUAL HMF'L I TUBE CASE 
GAMMA- 1 .38 00 OB 07 N= SO 
ITERATION SAMPLING PERIOD 

0 0 .GOO 00 00 GOOD- 01 

1 0.648 0879984 D- 01 

£ 0 . 6678086 0 1 9 It- 0 1 

3 0.67692321 2 OD- 01 

4 0.68 1 3 04 3 34 1 0 - 0 I 

5 0.63342504290-01 

6 0.63430395940-01 

7 0 . 684464554311- 0 1 

8 0 .6846072 173D- 01 


OF DISCRETE DESCRIBING FUNCTION 
K2®K3=0. 1 
PH I = 0.0 

AMPLITUDE 
0 . 1 00000000 OD- 04 
0.761 76c 9 8 54 D- 05 
0.717249318911-05 
0 .69968221 38 P- 05 
0 . 6 9 1 6*4 V J 1 5 0 D — 05 
0.6832531 139 D- 05 
0 • 6868 l 9 06 1 2D— 05 
0 . 6858 1 629 08 B- 05 
0. 685474493 OD- 05 


ND= 1 SAMPLING F‘ERIOD= 6 .84679D- 02. SEC AMPL I TUBE= 6.S5169D-06 N= 3 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1 .38 00 OD 07 N= 18 
ITERATION SAMPLING PERIOD 

0 O.6346737904P-O1 

1 0 . ? 043256 7 0 7 D- 0 1 

£ 0, 7274274 026D- 01 

3 0.7391 2561 5 ID- 01 

4 0 .7449677609D- 0 1 

5 0 .7478967216D- 01 

6 0.74 94 08586 IP- 01 

7 0 . 75 00 2 55 1 96D- 0 1 


□F DISCRETE DESCRIBING FUNCTION 
K2=K>0.1 

PHI= 0.0 

AMPLITUDE 
0 ■ 6*85 1 69266*2 D — U5 
0 . 7494237996D— 05 
0 . 73585732 78D- 05 
0 . 727 08 07448D- 05 
0 . 722 4 6 £ 022 2 D - 0 5 
0 .720121 953 9 D - 0 5 
0.71 3 92 69 556 D- 05 
0.7187 9 3 6 3 5 5 B - 0 5 


NB= 1 SAMPLING PERIOD® 7.5 0 03 ID- 02 SE C AMPLITUDE* 7.181 46 1 <- 06 N= 7 


LST SYSTEM -NUMER I ChL SOLUT I ON 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA = 1 .380 OOP 07 N= 16 
ITERATION SAMPLING PERIOD 

0 0. 75 0 031 3876 D- 01 

1 0.775 72 07225 D- 0 1 

2 0.8 0 3 6 5 2 6 5 2 3 D - 0 1 

3 0,8177101 684 D- 0 1 

4 0 . 824 7 1 75785D- 0 1 

5 0 . 8 28 22 0 38 5 7 U - 0 1 

6 0.8 3 0 0 1817 9 3 P - 0 l 

7 0.83091 001 45P- 01 

S 0 .83 092497 Q7D- 0 1 


OF DISCRETE. DESCRIBING FUNCTION 
K2=K3=0 . 1 
F'HI = O.O. 

AMPLITUDE 
0.718I455222D-05 
0 .7904 l 02565P - 05 
0 . 7 ? § 3 5 5 6 8 3 2 D -05 
0 . 76993299 1 2D- 05 
0 .7654 03 1462D- 05 
0 . 763 1 0 0342 1 P- 05 
0.7619 0625 86 D- 05 
0 . 76 1 567974 1 D- 05 
0 .761 0629 3 0 2 D ~ 0 5 


NP= 1 SAMPLING PERIOD®* 3 . 3 095 0 D - 0 £ SEC AMPLITUDE® 7.61094D-06 N® 3 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE K2=K3=0.1 


GAMMA® I .33 

0.0 OD 07 N® 14 

PHI® O.O 

ITERATION 

. SAMPLING PERIOD 

AMPLITUDE 

0 

0 . 83 095 0 1 073D- 0 1 

0 . 76 1 09376 09P- 05 

1 

i'i . 067 1 839 029 D- 0 1 

0 . 2:4 0364255 1 D- 05 

•"» 

L_ 

0 .901 7647994 D- 0 1 

0.331839 1 1 06 D- 05 

3 

0.9190691 383 D - 0 1 

0 . 8 2 4 6 3 0 7 4 9 8 D - 0 5 

4 

0.9276881881 D-01 

0 . 82 0 5510 5 5 9 D - 05 

5 

0.93 1986 12 04 P-01 

0.31 2-45 58899D- 05 

6 

0 . 934 1 75 0696 D- 0 1 

0.81 73583624 D- 05 

f 

0 . 935 34989 4 2 D - 0 1 

0 . 3 1 6 3 723 8 4 9 D - 05 

8 

0. 9353941 6 07D- 01 

0 . 8 1 67792323D- 05 


NP= 1 SAMPLING PERIOD® 9 .35426D-G2SEC AMPLITUDE® S.166S9B-06 N® 8 


Figure 2-7h 
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LST SYSTEM-NUMERICAL SOLUTION 
TUG AXIS-EQUAL AMPLITUDE CASE 
GAMMA- 1.38 00 OP 07 H= . IS 
ITERATION SAMPLING PERIOD 

0 0 .9354264774D-01 

1 0.99 ft. 1 IS888D-01 

2 0 . 1 03 5 4 39 72 3D 0 0 

3 (i . 1 OS' 7 S.'S “1 S6 SD 0 0 

-1 0-1 (I 3 74 79 [I 0 If 

5 0-1 073668 1 1 9 M ij 0 

A 0.10764 09 93 5 D 0 0 

7 0.1 077889937D i.iij 

. 8 0.1 078 16 0584 D 00 


OF DISCRETE DESCRIBING FUNCTION 

K2=K3=0.1 

PHI = 0.0 

AM PL I TUBE 
0.81668857990-05 
0 . 6 9 901568550-05 
0 .8984 l 68 0860- 05 
0 .8944 l 1987611- 05 
0 .8 -* I 6 74 *>9 (.5 0 • n5 
0 . 89 088 0 7 u 84 0- 0*5 
0 . 88 9 451 8 61 6 D - 05 
0 .8890091 778D- 05 
0 . 88 9 4 5 3 9 0 1 6 D - 0 5 


ND= 1 SAMPLING PERIOD* 1 . 07818D-01SEC AMPLITUDE* S.88944D-06 N= 8 


LST SYSTEM-NUMERICAL SOLUTION OP DISCRETE DESCRIBING FUNCTION 


TWO AXIS- 

-EQUAL AMPLITUDE t 

CASE 

K2=K3=0. 1 

GAMMA* 1 . 

38 0000 07 N= 

1 0 

PHI* 0.0 

ITERATION SAMPLING PERIOD 

AMPLITUDE 

0 

0.1 0781781 13D 

0 o 

0 . 88894 44839 D- 05 

1 

0.11771 00597D 

0 0 

0 . 962 141 54 77D- 05 

-i. 

0.12323362191' 

0 0 

0 .9301 3655 1 ?D- 05 

3 

0 . 186 05 038 4 3D 

0 0 

0 ■ 9 4 u t-9iL y y c D — o 5 

4 

0 . i 87 4 6 9 6 8 ijD 

0 0 

0 .934981 1926D— 05 

C 

■J 

0 . 128 167 033 OD 

00 

0 . 985224 5354 D- 05 

6 

0 . 1885 18898011 

0 0 

0 . 93536 1 1 699D- 05 

7 

0. 1887 044 7450 

0 0 

0 . 9 A 5 3 3 f - A * il 3D- 0 5 

8 

0 . 1838196333D 

0 0 

0 . 9 855 8 3 55 5 3D- 0 5 


ND= 1 SAMPLING PERIOD* 1 . 8882 0D- 0 1 SEC AMPLITUDE* 9 . 35437 D- 06 N= 8 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE K2=K3=0.1 
GAMMA* 1 .3800 011 n? N- 8 
ITERATION SAMPLING PERIOD 

0 0. 183S195875I) 0 0 

1 0.1 49737 1337D 00 

8 0. 15588 0S532D 00 

3 0 . 1593738457D 00 

4 0.16180071831! 00 

5 0.1681 1 35686 D 00 

6 0.1 685677 19 OD 00 

7 0.1688001 76 OH 0 0 

8 0.1 6894 753 OSD 0 0 

ND= 1 SAMPLING PERIOD* 1 .68946D- 0 1 SEC AMPLITUDE* 1 . 18823D- 05 N= 8 


PH I = 0 . 0 . 

AMPLITUDE 
0 .9854365774D- 05 
0.101 988965 0 D- 04 
0.1 08 04674 1 6D- 04 
0.11 08S85385D- 04 
0 . 1 1 186355 ?SD- 04 
0 . 1 1 1 7486 05 1 D- 04 
0.11 19856 75 7 D- 04 
0.1121 028633 D- 04 
0 . 1 1813291 15D-04 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 


TWO AXIS- 

EQUAL AMPLITUDE C 

RSE 

K2=K3=0. 1 

GAMMA* 1 . 

38 0 0 OD 07 N* 

ij 

PHI* 0.0 

I TERAT I DM S AMPL I MG PER I DD 

AMPLITUDE 

0 

0 • 1 6 8 9 4 w y I -7 c* 8 1' 

0 0 

0 . 1 1888894 l.l 9 Ii- 04 

1 

0.221 89675 06D 

0 0 

0 . 1 06597 1 96 OD- 04 

8 

0.8836 74 1 739H 

0 0 

0.11 98790559D- 04 

3 

0.28751 03138D 

0 0 

0. 1256581 187D- 04 

4 

0 . 2295949249 D 

0 0 

0. 1286281 146D-04 


0.23 065 735 I5D 

0 0 

0. 1301578755B-04 

6 

0.831 1396 028 D 

0 0 

0.13 093? 4889D — ij4 

7 

0.231 45705 19D 

0 0 

0. 131 33 13839 D- 04 

8 

0.23159677 8 3 D 

00 

0 . 1 3 1 5856779H- 04 

9 

0.231 71 93751 D 

00 

0 . 1 3 1 588 1 788B- 04 


ND* 1 SAMPLING PERIOD* 2 .3 1695D- 01SEC AMPLITUDE* 1.31578D-05 N* 9 


Figure 2-71 



LST SYSTEM-NUMERIC AL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE j K2= l.uGOD-Ol, K 3= 1.000D-01 


G AM M A = 1.33 0 0 0 D 0 7 N = 1 9 

PH I = 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0.69000 0 0 0 0 0 D - 0 1 

0 . 635 0 0 0 0 0 0 OP- 05 

1 

0 . 693 2 0 023 6 3 D - 0 1 

0.7! 043 79 736 D- 05 

2 

0.7067 193796D- 01 

0.7 06213372711-05 

3 

0.71 1 18131720-01 

0 - 7 0 323 06 3 0 0 D - 0 5 

4 

0. 7 1 3477636 I'D- 01 

0 .701 5467044D-05 

5 

0. 71 469586 03D- 01 

0.70 065 0 1 29 3D- 05 

6 

0.71 4962 0972 D- 01 

0.7 00774 11 75D- 05 


ND= 1 SAMPLING F'ERIOD= 7 . 14990D-02SEC AMPL I TUDE= 7.00315D-06 


LST SYSTEM- 

NUMERICAL SOLUTION 

OF DISCRETE DESCRIBING FUNCTION 

TWO AXIS-EQUAL AMPLITUDE CASE 

K2* 1 . 0 0 00-01 .» K3= 1 .0 000-01 

GAMMA* 1 .33 

00 OH 07 N= 17 

PHI* 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0.71 49897 4 c' 4 D— 0 1 

0.70031487640- 05 

1 

0.735962 122 ID- 01 

0 . 7 7 02 5 3 1 4 92 0 - 0 5 

■:> 

L. 

0.761 3679 043 D - 0 1 

0 ,75657762320-05 

3 

0.7 74 136 035 ID- 01 

0 .74 7791 01 740-05 

4 

0.78049 33 3 52 D - 0 1 

0 .7431759493D-05 


0.73367 1532 ID- 01 

0 . 74 084396820-05 

6 

0.7353 068334 D- 01 

0 .73965737340- 05 

7 

0.786 038523: 0 D -01 

0.73944889100-05 

NO* 1 -SAMPLING PERIOD* 7 .36051 D-02SEC AMPLITUDE* 7. 3332 ID- 

LST SYSTEM- 

NUMERICAL SOLUTION 

OF DISCRETE DESCRIBING FUNCTION 

TWO AXIS-EQ 

UAL AMPLITUDE CASE 

, K2* 1 . 0 0 0 D -01 * K 3= 1 .0 0 0 D-01 

GAMMA* 1.32 

fiOOD 07 N= 15 

PHI* 0.0 

ITERATION 

SAMPLING PERIOD 

HMPLITUDE 

0 

0 . 786 05 1 3935 D- 0 1 

0 . 73832 05 1 250- 05 

1 

0.8131 3597230- 0 1 

0 .3 1 91 99 0949D- 05 

O 

0 . 84 4 8 3 9 7 6 3 6 D - 0 1 

0 . 3 07095 086SD- 05 

3 

0 . 8593 065953D- 0 1 

0 .793725639 ID- 05 

4 

0 . 8 6 7 5 1 6 5 7 4 5 D - 0 1 

0 .7942 3 0 15070- . 0 5 

C 

0 . 37 1 3578826 D- 0 1 

0 .791 963 0257D- 05 

6 

0 .3 733 1 87 734 D- 0 1 

0 . 7907942364 0- 05 

7 

0 . 874 3 2:4 5 62 3 D- 0 1 

0.79038 5 6 5 5 3 P - 0 5 

c* 

0.3743 53 3233 D - 0 1 

0 .79002743 83 0 - 05 

ND= i SAMPLING PERIOD* 3 .74330B-02SEC AMPLITUDE* 7.90034D- 

LST SYSTEM- 

NUMERICAL SOLUTION 

OF DISCRETE DESCRIBING FUNCTION 

TWO AXIS-EQUAL AMPLITUDE CASE 

j K2= 1.0000-01* K3= 1.000D-01 

GAMMA* 1.35 

0 0 00, 07 N= 13 

PHI* 0.0 

ITERATION 

SAMPLING PERIOD 

AMPLITUDE 

0 

0 . 8 74 3 79 6 535 D~ 0 1 

0.790 0 3 3 8 7 02 D -05 

1 

0.91 175 154 780- 01 

0.832361 1 04 0D- 05 

2 

0 . 95 03524 084D- 0 1 

0 . 8739792 1 8 0D~ 05 

o 

0 . 970235 0537 D- 0 1 

0 . 8 6 6 3 9 9 0 2 3 0 P — 0 5 

4 

0 . 979859686 0D- 0 ! 

0 .36287 3 0 9 7 3 D - 0 5 

5 

0 . 98464 07971 D-01 

u • o 6 IJ I? 3 -J’ *4 *r* ill f D 0 

6 

0 . 987059 1 6420-0 1 

0 . 859 78 56 68 7 D- 05 , 

7 

0 .938 .>6189770- 0 1 

0 , 859276 1851 D* 05 

n 

■_> 

0.98344635570-01 

0.35942296 ISO- 05 


. NO* 1 SAMPLING PERIOD* 9 .83479D-02SEC HMPLITUDE* 3 .59222B-06 


Figure 2-7j 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE ♦ K3= 1 . 0 ODD- 01 > K3= 1 .000D-01 


GRMMA= 1 .38 00 OD 07 N= 11 PHI= 0.0 


ITERATION 

SAMPLING PERIOD 

amplitude 

0 

0 . 9b'84 ?94b'6bD- 

0 2 

0 .85933193788-05 

1 

0 . 1 0447474998 

00 

0 . 964 030 045 OD- 05 

3 

0. 109S 1373 798 

0 0 

0.9645574868D-05 

3 

0. 1 1801481 458 

00 

0.96093513548-05 

4 

0.11 325690738 

0 0 

0 . 9 58 3 7334 7 8 8 - 05 

5 

0.11 337364 1 9 D 

0 0 

0. 957 01 55354 D- 05 

6 

0.1 14131 14158 

0 0 

0 . 9 563 3853 048-05 

7 

0.1 14 3453988 D 

nn 

0 . 95593765 1 OD- 05 

S 

0 . 1 14 3983886 D 

00 

0 .9564537 l 08 D- 05 

9 

0. 1 1 4 399 1 88 7 D 

0 0 

0 . 9 554 8 7 933 0 D - 05 

10 

0.1 144 00643 OD 

00 

0 . 9 55 3 3 1 0 5 7 7 P - 0 5 


NP= 1 SAMPLING PERIDP= I . 144 0 ID- 01 SEC AMPL I TUDE= 9. 557438-06 


LST SYSTEM-NUMERICAL 


SOLUTION OF DISCRETE DESCRIBING FUNCTION 

1 . 00 OD- 01 


TWO AX I S- 

EQUAL AMPLITUDE CASE 

-• K3= 1 . 0 008-01 • K3 

GAMMA= 1 . 

38 0 0 OD 07 N= 

*5* 

PHI= 0.0 

ITERATION SAMPLING PERIOD 

AMPLITUDE 

0 

0.1 14401303 ID 

00 

0 .9557425397D- 05 

1 

0. 134 1494 01 7D 

00 

0 . 1 0677543938- 04 

3 

0.13065877698 

0 0 

0 .1 091 30 1 3358-04 

3 

0 . 1 3395 0732 OD 

00 

0 . 1 09742 09568- 04 

4 

0 . 135604 0143D 

00 

0.1 099 0449338-04 

5 

0.1 364338 1338 

00 

0 . 1 09974679 i p- 04 

6 

0.136331 14 16D 

00 

0,1 1 00 1 3556 1 D- 04 

-7 

0 .13704193198 

00 

0.11 0 0364 376 8- 04 

0 

0. 137 1781 S79D 

00 

0.11 003694 04 D- 04 


ND= 1 


SAMPLING PER I OD= 1 .371738-01 SEC 


AMPL I TUDE= 1.1004£D-0: 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 


TWD AXIS-EQUAL AMPLITUDE CASE 
GAMMA = 1 .380008 07 N= ' 7 
ITERATION SAMPLING PER I OH 


K2= 1 . OOOB-Ol 
PHI = 0.0 

AMPLITUDE 


K'3= 1. OOOD-Ol 


0 

0 . 1 37 1 78 1 S3 3D 

0 0 

0.1 1 0 04 1 99398- 04 

1 

0 . 1 533943 136D 

00 

0.11 836904388- 04 

3 

0 . 1 66 1 3596368 

0 0 

0.1 

8743347538-04 

3 

0.1704 1S3536D 

00 

0.1 

3 07466534 D- 04 

4 

0.1 7363638368 

0 0 

0.1 

3 2 3 0 5 10068-04 

5 

0,17374430758 

00 

0 . 1 

3293005778-04 

. 6 

0.1 74393797 ID 

on 

0.1 

333879 1 868- 04 

7 

0.17456944338 

0 0 

■0.1 

3347131518-04 

8 

0.17473399338 

0 0 

0 . 1 

3 3 5 4 08 0 5 6 8 - 0 4 

9 

0.1 74558 13378 

0 0 

0.1 

3370738758- 04 

1 0 

0.17474059658 

nn 

0. 1 

3 36 9 2 6 5 3 5 8-04 

1 1 

0.17474079518 

0 0 

0.1 

3363517338-04 

NP= 1 SAMPLING PER IDD= 1. 

7474 

3D- 01 

SEC AMPLITUDE^ 

LST SYSTEM- 

NUMERICAL SOLUTION 

OF DI 

SCRETE BESCRIBINt 

TWO AX IS- EQUAL AMPLITUDE C 

ATE , 

k 8 = 

1 . uOOD-Ol > K3= 1 

6AMMA= 1 .33 
ITERATION 

000D 07 N= 5 

SAMPLING PERIOD 

PH I = 

0 . 0 

amplitude 

0 

0. 17474333658 

0 0 

0.1 

2 26 8 1 7 5 5 5 D — 0 4 

1 

0 .24 9 5*2 96t' 5b P 

00 

0 . 1 

875924454 P- 0 4 

O 

0 .24 6 ( 0692 4 9 D 

00 

0. 1494333451 D- 04 

3 

0.351 03 193 058 

0 0 

0.1 

59 5 8 38 8 7 3 D - 0 4 

4 

0 .3535531 7358 

0 0 

0.1 

6497478 1 9 D- 04 

5 

0.35484561648 

0 0 

0.1 

6- 7 8 31 6481 8— 0 4 

6 

0.35549593188 

00 

0 . 16930855498-04 

7 

0. 8558886254 D 

0 0 

0.1 

? 005 00828 D- 0 4 

0 

0.85593943168 

0 0 

0.1 

7 04 85 0633 D- 04 

9 

0 .85608816688 

0 0 

0.1 

7 06 0 1 5366 8- 04 

10 

0 .85603333848 

0 0 

0.17 07718967D-04 


ND= I SAMPLING PEPIDD= 3 .560848-01 SEC AMPLITUDE^ 1 .70779D-0! 


N= 1 0 


N= 3 


N= 1 1 


N= 1 0 


Figure 2- 7k 
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LST SYSTEM-NUMERICAL SOLUTION 
TUD AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .38 GOOD 07 N= 4 


ITERATION 

SAMPLING PERIOD 

0 

0 .5000 00 000 OD 

00 

1 

0.5 15424 0480D 

0 0 

2 

0 . 5395948 07 OD 

00 

3 

0 .548 0038 734 0 

0 0 

4 

0.5501 2329360 

0 0 

t 

0.54983533 OSD 

0 0 

6 

0.54961532560 

0 0 

i 

0.54955191 1 ID 

0 0 

8 

0.54955726570 

00 


OF DISCRETE DESCRIBING FUNCTION 
K2= 1 .OGOD-Ol ’ K3= 1.000H-01 
PH I = 9.0 000 OD 01 
AMPLITUDE 
0.15000000000-04 
0. 137790436 1 D- 04 
0.14465672950-04 
0 . 1 4 7 4 7 8 1 5 0 9 D - 04 
0 . 1 4869462980- 04 
0 . 1 4933778980- 04 
0 . 1 4 96763 755 B- 04 
0.14985359140-04 
0 . 14981260900-04 


NO* 1 SAMPLING PERIOD* 5 .496370- 01 SEC AMPLI TUBE= 1 .499130-05 N= 8 


LST SYSTEM-NUMERICAL SOLUTION 
TWO AXIS-EQUAL AMPLITUDE CASE 
GAMMA* 1 .380000 07 M- 4 
ITERATION SAMPLING PERIOD 

0 0 .54968655340 00 

2 5.496370-01 -1.91 2740- 

2 5.5 01 360- 01 -1.956880- 

1 0.4 234 1 86 7960 0 0 

2 4 . 234 1 90-0 1 -7. 38 1 8 OD- 

£ 4.238420-01 -7.458110 

2 0.-37 3 0 698 8 9 1 0 0 0 

3 0.41284341440 00 

4 0.431 8546 1 790 • 0 0 

5 0.43974015240 00 

6 0.44 3 83 0 1096 D 0 0 

7 0.44584 33 1951) 0 0 

S 0 .44685247930 00 

9 0 .44739 068880 0 0 

10 0.44752056080 00 


OF DISCRETE DESCRIBING FUNCTION 
K2= 1 . 0 0 OD- 0 1 , K3= 1 . 0 0 OD- 0 1 
PHI* S.00000D 01 
AMPLITUDE 
0. 14991274 780-04 
05 
05 

0 . 845 4 5633 0 3 D - 0 4 
05 
- 05 

0 .16924221450-04 
0 . 1 6376256940-04 
0 .1 6549076560-04 
0 . 16651 229 190- 04 
0. 16-704 043940- 04 
0 . 1 673 0648470- 04 
0 .1674 39 08330- 04 
0.16750107190-04 
0.1 675 157 0260- 04 


NO* 1 SAMPLING PERIOD* 4. 47512D-01SEC AMPLITUOE= 1.674880-05 N=10 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE , K2= 1.0000-01 , K3= 1 .0000-01 


GfthMA= 1 .38 

0-000 07 N= 

4 

PHI = 

7.000000 01 

ITERATION 

SAMPLING PEP IDO 


AMPLITUDE 

0 

0 .44751201780 

0 0 

0 . 

167488 t 54 00-04 

1 

0.42160824570 

0 0 

o . 

186 1 1 704290- 04 

c! 

0 .405 1 0479020 

00 

o . 

1 '? 0 0 0 o c 1 *4 4 B 04 

3 

0.39333940790 

oo 

0 . 

19183136370-04 

4 

0 .3367604 1540 

0 0 

0 . 

19281807400-04 

5 

0.38254075020 

00 

0 . 

1 9 3352545 9 D - 0 4 

6 

0.3301 1 777640 

0 0 

0 . 

19263364770-04 

? 

0. 3787622 054D 

0 0 

0 . 

19377729050-04 

S 

0.37801602880 

0 0 

0 . 

19384922950-04 

o 

0. 377607655 OD 

0 0 

- 0 . 

1938848 19 OD- 04 

1 0 

0.37738238440 

00 

0 . 

19390234 1 3D- 04 


ND= 1 SAMPLING PERIOD* 3 .77255D-01SEC AMPLITUDE* 1. 939110-05 N=10 


Figure 2 -8a 
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LST SYSTEM-NUMERICAL SOLUT I ON OF HIS CRETE DESCRIBING FUNCTinN 
TWO AXIS-EQUAL AMPLITUDE CASE? K2= 1.000D-01* K3= 1 . 0000-01 
GAMMA® 1 .33 GOOD 07 N = 4- F'H l = 6.00000D 01 

ITERAT ION SAMPLING PERIOD AMPLITUDE 

0 0. 37725499420 00 0 . 1 9391 1 0275D-04 

1 O .3599882449D no 0 .212231 3759D- 04 

2 0 .34 ?4431 1 5 l D 00 0 .2 1994930760- 04 

3 0.3339745293D 00 0 .224639471 3D- 04 

4 0.3336445753B 0 0 0 .227517122411-04 

5 0.3304527295D 0 0 0.22921154360-04 

6 0 .3286099425B 00 . 0 .230 1675458D-04 

? 0 .32757344*. ID 0 0 0.23063949230-04 

8 . 0.327000 1 735B 00 0 .2309677767D-04 

9 0 .32663529320 00 0 .231 1 134342D-04 

NB= 1 SAMPLING PERIOD® 3 .265 11 D- 01 SEC AMPLITUDE® S.31189B-05 N= 9 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AX I S-EQUAL AMPLITUDE CASE* K2= l.OQOD-Ol* K3= 1. GOOD- 01 
GAMMA® 1. 380 0 OIi 07 N= 4 F‘HI= 5.0 00 0 0D 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 . 3265 1 1 468 OH 0 0 0 . 23 1 1 S85352D- 04 

1 0.313443703811 00 0 .24969591 09D-04 

2 0.304 3395 09 5 D 00 0 .2615 63 8 5 1 0 D - 0 4 

3 0 ,293 36 02 071 D 0 0 0 .2692 65383 6 D - 0 4 

4 0 .2946473529D 00 0 .274 055606 OB- 04 

5 0.29244 0829 3 D 0 0 0 . 276865 9 223 D - 0 4 

6 0 .291 171 7454B 00 0 .2734892 1 57B- 04 

7 0 .2904579278D 00 0 .27937122550- 04 

8 0 .290 06 0381 2B 00 0 .279849304 0D- 04 

9 • 0 . 289836 G556B 0 0 0.2301 1 33966D- 04 

ND= 1 SAMPLING PERIOD® 2 .89699D-01SEC AMPLITUDE® 2.802850-05 N~ 9 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE C ASE * K2= 1 . GOOD- 01 » K3- 1.0 0 OB- 01 
GAMMA® 1.38000B 07 N= 4 PHI® 4.000000 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0.289699014 9 D 00 0,28028478030-04 

1 0 .27948906780 00 0 .300 1 1331 030- 04 . 

2 0. 27299 17247D 00 0.314641 03 14D- 04 

3 0.26 8 9976 36 9 B 00 0 .38426271840-04 

4 0.2 6 6- 623 51 140 00 0 . 3 3 0 1 5 1 3 4 6 1 0 - 04 

5 0 , lit- 5* 258468 30 UO 0 . 332557499GB— U4 

6 0.26447724300 00 0.33545739330-04 

7 0.26404811820 00 0 . 33650275380-04 

8 0.263 78963 53 0 0 0 0.3371 1 0854 9B- 04 

9 0 . 26 35933842 0 0 0 0 . 3 37714 4 3 660- 04 

10 0 .26353728 1 OO 00 0 .33784295770-04 . 

11 0.26346992600 00 0.33823587680-04 

12 0 . c.*68 36*42669 0 0 0 m r 4 5 2 o 7520 — u4 

13 0 . 36 3 3 1 0 1 3 9 9 1* 0 0 0 . 3 39 4466571B - 0 4 

14 0 .263573 0871 D 0 0 0.337 04942450-04 

15 0.26357300690 00 0 .336804 0 1 7 1 0-04 

ND= 1 SAMPLING PERIOD® 2 . 635730-0 1 SEC 'AMPLITUDE® 3.36804D-05 N=15 


Figure 2-8b 
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LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* E£= 1. 0000-01, K3* l.GCinD-01 
GAMMA* 1.38 00 OP 0? N= 4 PHI* 3. 000 00D Ut 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0. £63573006?® 00 0.3363 04 01 7 ID- 04 

1 0. 25572 0G863D 00 0 .3568292423P-04 

2 0.2512S15177B 00 0 .3709936920D-04 

3 0.248801 13S7D 00 0 . 37995265 l 3D- 04 

4 0.2474297532D 00 0.38519499431-04 

5 0 .24 66 7 7324 3D 00 0.3881 19P135D-04 

6 .0. 24626234 Q9B 0 0 0 .389725621 7D— 04 

7 0 . £46 0 1 7296 CD 0 0 0 . 39 07037428D- 04 

S 0 .24561 14935D 00 0 .39348503661'- 04 

9 0. £4597 1S36 Oil 00 0 .3364238696D-04 

10 0 .£45831 0873P 00 0 .389 165933 0D- 04 

11 0 . £4 5 74 338 56 It 0 0 0 . 3 9 04 8 0 1 4 34 D- 04 

12 0 .£46 0 09327 00 00 0 .39 11 03742 IB-04 

13 0.. £457639879D 00 0 . 3 92 4 2 92 7 £7 D- 04 

14 0 .£4533068470 00 ' O .398189791 ID-04 

15 0 . £45643 042 3 Ii 0 0 0 . 393 Q859796D— 04 

MD- 1 SAMPLING PERIOD* £ .45643D-01SEC AMPLITUDE* 3.930S6D-05 M*15 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AX I S- EQUAL AMPLITUDE CASE > K£= 1.000D-01* K3= 1.000D-01 
■GAMMA* 1 .38 0 0 OD 07 N= 4 PH I - £.000001* 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 .245643 0421 Ii 00 0 .3930859796D-04 

1 0 .24034 1 5437D 00 0 .40393445 0£H- 04 

£ 0 . £sy £26*2^* 1 4 D 00 fj ■ 4 U t. S' 6 c* c o .it o D — U4 

3 0.2374380 085B 00 0 .4064485320D-04 

4 0.2371 31 5275D 0 0 0 .4 0549597 08D- 04 

5 0 . £36978 1 829B 0 0 0.4 05 0 046427D- 04 

6 0. £368591 6 091* 0 0 0.4 052 03689711-04 

ND= 1 SAMPLING PERIOD* £ .36794 D- 01 SEC AMPLITUDE* 4 . 054 1 2D- 05 N= 6 


LST SYSTEM-NUMERICAL SOLUTION OF DISCRETE DESCRIBING FUNCTION 
TWO AXIS-EQUAL AMPLITUDE CASE* h 2= l.OOOD-Ol* K 3- I.000D-01 
GAMMA* l .380 00 D 07 N* 4 PHI= 1.00 00 OD 01 
ITERATION SAMPLING PERIOD AMPLITUDE 

0 0 .322 3 7 9 9 3 5 8 Ii 0 0 0 . 167469409 3 D ~ 0 4 

1 0 . 3 8 1 7 32' 7 9 8 5 D 0 0 0.8144 6 0 2 4 1 6 1 ' - 0"4 

2 0 . 3 479304 1 3 3 D 0 0 0 . 1476 8 5 1 8 £60-04 

3 0 - 3 5 6- 0 0 4 4 3 9 0 D 0 0 0 . 146 9 513 7 0 6 D “04 

4 0 . 3 5 9 3 4 3 0 0 0 3 0 0 0 0 .14684 8 7 7 £ 8 D - 0 4 

MP* 1 SAMPLING PERIOD* 3 .60729H-01 SEC AMPLITUDE* 1.46858D-05 NT* 4 


Figure 2-8c 





(sec) 


34 



Amplitude ( * 10”^) 

Figure 2-10. Amplitude and frequency (2ir/nT) of self-sustained oscillations 
for various sampling periods in the discrete-data two-axis LST 
system; y - 1.38xl0 7 . 



(sec) 
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Amplitude ( * 10 - ' 7 ) 

Figure 2-11, Amplitude and frequency ( 27 t/ nT ) of self-sustained oscillations 
for various sampling periods in the discrete-data two-axis LST 
system; y = 1 . 38*1 O' 7 - 
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